Spin-Peierls phases in pyrochlore antiferromagnets 
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In the highly frustrated pyrochlore magnet spins form a lattice of corner- sharing tetrahedra. 
We show that the tetrahedral "molecule" at the heart of this structure undergoes a Jahn- Teller 
distortion when lattice motion is coupled to the antiferromagnetism. We extend this analysis to the 
full pyrochlore lattice by means of Landau theory and argue that it should exhibit "spin-Peierls" 
phases with bond order but no spin order. We find a range of Neel phases, with collinear, coplanar 
and noncoplanar order. While collinear Neel phases are easiest to generate microscopically, we also 
exhibit an interaction that gives rise to a coplanar state instead. 

PACS numbers: 75.10.Hk, 75.10.Jm. 



I. INTRODUCTION 

The study of highly frustrated magnets began with 
Wannier and Houtappel's realization that the triangu- 
lar lattice Ising antiferromagnet is paramagnetic at any 
nonzero temperature and exhibits a macroscopic entropy 
even in at zero temperature.EII3 This canonical exam- 
ple illustrates the defining characteristics of such sys- 
tems - their failure to order at temperatures of order 
the exchange constant, empirically derivable from the 
high temperature Curie susceptibility, and a large low- 
temperature entropy.Erea 

The advent of the cuprate superconductors led to seri- 
ously renewed interest in these systems in the hope of 
finding a quantum spin liquid - the zero temperature 
state of a quantum magnet that fails to order. Subse- 
quently their study has blossomed, driven by an increas- 
ing list of materials that exhibit highly frustrated antifer- 
romagnetism, and is driven as much in hopes of finding 
unusual ordering at low temperatures. An appealing, if 
optimistic, analogy is to the quantum Hall system, where 
the magnetic field frustrates the kinetic energy and pro- 
duces a macroscopic degeneracy, which is then lifted by 
residual terms in the Hamiltonian to produce a rich phase 
diagram with various orderings. 

The most promising system in this regard is the 
nearest-neighbor Heisenberg system on the "pyrochlore" 
lattice, a network of corner-sharing tetrahedara (Fig [I]). 
The idealized system has a vast ground state degener- 
acy in the classical limit of infinite spin and there is 
a large and growing list of materials that approximate 
this to varying degrees, including doped variants that su- 
perconduct or display behavior reminiscent of the heavy 
fermions. 

In this paper we limit ourselves to the insulating mag- 
nets. Here it is knownB that the classical system does not 




FIG. 1: The "pyrochlore" lattice. Magnetic ions are situated 
in the corners of tetrahedra. 



exhibit order by disorder and remains a (co-operative) 
paramagnet down to T — 0. On general grounds one ex- 
pects that quantum fluctuations will select an ordering 
at sufficiently large spin in a, spin- wave treatment about 
the classical ground states.H At smaller values of spin, 
the situation is unsettled with some form of singlet order 
likely-Ha 

Recently, building on work of Yamashita and Ueda,tJ 
we have shown that in the presence of a coupling to 
the lattice, a different mechanism of degeneracy lifting is 
likely to operate.t£l This involves a version of the Jahn- 
Teller effect in which the lattice distorts to gain exchange 
energy (the "spin- Teller" effect) and thereby relieves the 
frustration. As magnetoelastic couplings are ubiquitous, 
and lead to a transition even at infinite spin, this mech- 
anism will dominate over the purely quantum selection 
effect, likely starting at modest values of the spin.EJ Also 
noteworthy in this problem is the likelihood of a finite 
temperature bond ordered phase preceeding the eventual 
establishment of Neel order. 

In this paper we present a detailed account of our ana- 
ysis of the Jahn- Teller effect for Heisenberg magnets on 
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the pyrochlore lattice. Parts of this work have already 
been summarised in a short paper .Ej In Section || we be- 
gin with the symmetry analysis of the Jahn- Teller distor- 
tion of a single tetrahedron in the classical limit and then 
extend it to q = phonons for the infinite lattice. Hav- 
ing identified the order parameter for lattice distortions 
(bond ordering) in this fashion, in Section III we con- 



struct a Landau theory of the transition into the bond or- 
dered state which we contrast with the Landau theory of 
the spin-Peierls transition in quasi-one dimensional sys- 
tems. Finally we turn to the establishment of Neel order 
which we discuss with the insight gained from analyzing 
bond order (Section [V). Such ordering is most natu- 
rally collinear but symmetry permits coplanarity and in 
Appendix ^ we describe an interaction that would bring 
it about. Appendix [b| gives the quantum theory of the 
Jahn- Teller effect in a single tetrahedron. The experi- 
mental situation with regard to structure is briefly dis- 
cussed in our concluding remarks in Section |v|. We will 
discuss the dynamical signatures of various phases in a 
forthcoming publication. 
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FIG. 2: Four spins of a tetrahedron in a ground state (zero 
total spin). 



II. JAHN-TELLER EFFECT 
A. Single tetrahedron 

The structural unit of the pyrochlore antiferromagnet 
is a tetrahedral "molecule" with four spins in the corners. 
Its high symmetry and the degeneracy of the ground state 
are the two prerequisites for the Jahn- Teller effect: the 
tetrahedron is distorted in its ground state. The ten- 
dency of individual tetrahedra to deform induces a co- 
herent distortion of the entire crystal. We will describe 
the Jahn- Teller effect for a single tetrahedron in detail to 
understand which aspects are relevant for the description 
of the spin-Peierls effect on the entire lattice. 

The energy of four spins on a regular tetrahedron is 

E = JJ2 s * ' s j ■ = o (Si + S 2 + S 3 + S 4 ) 2 - 2 JS(S + 1). 

i<j 

(1) 

In a ground state, the total spin is 0. For quantum spins 
of length S, there are 25 + 1 linearly independent ground 
states, which can be constructed as follows. The total 
spin of the pair Si and S 2 can be 0, 1, 2, . . . , 2S, and like- 
wise for the other pair, S3 and S4. An overall spin singlet 
can be formed by combining two singlets, two triplets, 
and so on, giving a total of 2S + 1 physically different 
singlet states. 

The problem of quantum spins on an elastic tetrahe- 
dron can be solved straightforwardly and is treated in 
detail in Appendix |B[ As the degeneracy and hence the 
Jahn- Teller distortion survive at arbitrarily large values 
of spin, the outcome (with the exception of the extreme 
quantum case S = 1/2) can be understood by looking at 
the simpler problem with classical spins. In essence this 



spin-Peierls effect is classical, in contrast with the usual 
cases where it goes away in that limit. 

For classical spins (S — > 00), the degeneracy of the 
ground state becomes continuous. In addition to a trivial 
rigid rotation of all four spins, there are two parameters 
that can be used for characterization of a ground state 
(Fig. |^) : the angle 28 between spins 1 and 2 and the angle 
2(f> between the planes 12 and 34. These two parameters 
determine the bond variables in a ground state, 

Si ■ S 2 = S 3 ■ S 4 = S ,2 cos26», 

S 3 - 83 = 81-84 = S 2 (sin 2 6 cos 20- cos 2 9), (2) 
S 3 -Si = S 2 -S 4 = — S 2 (sin 2 6 cos 2(f) -(- cos 2 9) . 

At the heart of the effect lies the dependence of ex- 
change interaction on the relative positions of spins. For 
example, if depends strictly on the distance between 
spins i and j, the contribution of this pair to exchange 
energy, 

Ekj — [J + {dJ/dr)5 nj + . . .](S; • Sj), 

generally has a term linear in the relative displacement 
Srij. Therefore, spins i and j exert on each other a force 
— (dJ/dr)(Si-Sj), which is repulsive or attractive depend- 
ing on the angle between the spins. In a generic ground 
state (Fig. ||), angles between spins are unequal, so that 
disparate forces cause a deformation of the tetrahedron. 

More generally, exchange interaction may depend not 
only on the distances between spins, but also on the an- 
gles between the bonds. We therefore write the magnetic 
and elastic energies of the spins in the most general form: 

E = E + J~](dJij/dx a )(Si ■ Sj)x a + ^2 kabXaXb/2- 

a,i,j a, b 

(3) 
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FIG. 3: Vibrational modes of a tetrahedral molecule. 

Here Eq is the energy of a ground state, x\ . . . X\<x are 
Cartesian coordinates of the spins and fc t, are the ap- 
propriate elastic constants. To reduce the number of 
independent coordinates and forces, it is convenient to 
classify them in terms of irreducible representations of 
the symmetry group Td of the tetrahedron. Using the 
appropriate linear combinations of coordinates x pa and 
forces f pa (where p labels irreducible representations, a 
enumerates its components) one obtains a simpler result: 

E = E Q +Y J [-J' P f Pa x pa + k p x 2 pa /2]. (4) 

Six vibrational modes may affect the exchange energy: 
a singlet A\, a doublet E, and a triplet T 2 . The breath- 
ing mode A\ uniformly rescales exchange interactions on 
all bonds and does not discriminate between different 
ground states; therefore it can be left out of consider- 
ation. A component of the vector triplet T 2 stretches 
and contracts by the same amount two bonds opposite 
each other (Fig. |J). As can be inferred from Eq. (||), 
such bonds are equally satisfied (or equally frustrated) 
in any ground state. Therefore stretching one and con- 
tracting the other to the same extent cancels the linear 
term in magnetic energy, making the triplet mode inef- 
fectual in relieving frustration via the Jahn- Teller mech- 
anism. Only one irreducible representation causes the 
Jahn- Teller effect: the doublet E whose components are 
tetragonal and orthorombic distortions of the tetrahe- 
dron (Fig. ||). Since no other representations will be dealt 
with, we will suppress the representation subscript p = E 
in what follows. 

The six bond variables S^ • Sj contain the same repre- 
sentations. The singlet A\ is the symmetric sum, which is 
nothing but the energy of the undistorted ground state 
([j]), i.e., a constant that does not favor any particular 
ground state. The triplet T 2 contains the differences of 
forces on opposite bonds, 

Si • S3 — S 2 • S4, Si • S4 — S2 • S3, Si • S2 — S3 • S4. 

(5) 




FIG. 4: The domain of the bond vector f = (/lj/z) = 
(/cos#, / sin#) is bounded by an equilateral triangle in the 
(fijfa) plane. Also shown are six extremal spin configura- 
tions. Strong (weak) bonds are denoted by solid (dashed) 
lines. The color of a state is determined by the color of frus- 
trated bonds. 

As already mentioned, these differences vanish in a 
ground state. The remaining forces form a doublet E 
showing the disparities between adjacent bonds: 

ft = [(Si + S 2 ) • (S 3 + S 4 ) - 2Si ■ S 2 - 2S 3 • s 4 ] /Vu, 
h = (Si - S a ) • (S 3 - S 4 )/2. (6) 

The component fx shows by how much the blue bonds are 
stronger than the rest (Fig. [IJ) ; / 2 measures the difference 
between the red and green bonds. The domain of possi- 
ble values of the vector f = / 2 ) = (/ cos a, / sin a) is 
an equilateral triangle. Its perimeter is made of coplanar 
ground states; the three corners correspond to the three 
distinct collinear ground states; they are marked by the 
color of frustrated bonds. The two components of f , like 
the two angles in Fig. |^, can be used to parametrize de- 
generate classical ground states of a tetrahedron; in fact, 
/1 = 2S 2 (1 - 3 cos 2 0)/y/3, h = 2S 2 sin 2 6 cos 2c/.'. 

After these simplifications, the energy of the system 
has the form 

E = E - J'f-x + fc.x 2 /2, (7) 

where x = (xi,x 2 ) are amplitudes of the tetragonal and 
orthorombic distortions. J' and k are the magnetic and 
elastic constants appropriate for the E representation. 
The energy is minimized when /ex = J'f , so that 

£min = E - J' 2 f/2k. (8) 
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FIG. 5: The Neel state obtained in the magnetoelastic model 
with q = phonons. The E g phonon mode dominates. The 
distortion weakens blue bonds on all tetrahedra (color scheme 
of Fig.g). 

One can view the — f 2 term as a quartic spin interaction 

~^T = ~C ^ (Si ' Sj)2+ const (9) 

induced by "integrating out" the phonons .0 It evidently 
prefers collinear ground states with four maximally sat- 
isfied bonds and two maximally frustrated bonds. The 
resulting distortion of the tetrahedron is tetragonal. It 
flattens or elongates along one of its C2 axes, depending 
on the sign of the derivative J'. 

Modulo global rotations of the spins, there are three 
degenerate ground states, which we label with the colors 
red, blue and green according to which pair of opposite 
bonds is frustrated (Fig. Their opposites are cyan, 
magenta, and yellow states with four frustrated bonds. 
We have found that ground states can have these sec- 
ondary colors in a model with more general spin interac- 
tions, e.g., 4-spin cyclic exchanges. See Appendix |a| for 
details. 

A consideration of quantum spins in the Appendix [b| 
yields essentially the same result: the energy is mini- 
mized when two opposite bonds (e.g., 12 and 34) have 
the highest spins 25* each. Such states are the quantum 
analogue of parallel spins. In contrast to the spin-Peierls 
effect on a Heisenberg chain, this one is a classical af- 
fair: instead of forming spin singlet on stronger bonds, 
Heisenberg spins of a tetrahedron form the highest spin 
on two weak bonds. 




FIG. 6: Same as Fig. q but with the dominant E u phonon. 
Red and green bonds are weakened alternatively. 



B. Pyrochlore lattice: q = phonons 

An attempt to extend this calculation to an infinite 
network of tetrahedra runs into a substantial problem: 
all possible phonon modes, the number of which is pro- 
portional to the number of tetrahedra, couple to bond 
variables, and any one of them may thus lead to a magne- 
toelastic distortion. In order to describe the basic physics 
of the Jahn- Teller effect on the full lattice, in this sec- 
tion we restrict ourselves to the (over)simplified version 
with only a few participating phonons. As a result of 
this crude approximation, "integrating out" the phonons 
produces an infinite-range interaction between vectors f 
of different tetrahedra. In a realistic model including 
phonons of all wavelengths, such forces will have a finite 
radius. However, the structure of the ground state is 
often insensitive to such details. 

We specialize to the case of phonons with the wavevec- 
tor q = 0. In effect, this restricts all tetrahedra of the 
same orientation to have the same distortion (Figs. and 
||) . The existence of two types of tetrahedra (labeled in 
what follows A and B), which differ by orientation, is 
the only new degree of freedom. The symmetry group 
of the lattice (with equivalent tetrahedra identified) is 
extended from Td to I x Td = Oh by the operation of in- 
version on any site, which exchanges tetrahedra A and B. 
Irreducible representations are those of Td labeled by an 
additional quantum number, parity under the inversion. 
The relevant phonons are E g and E u , which are, respec- 
tively, uniform or staggered distortions of the lattice. For 
example, the first component of E g stretches all tetrahe- 
dra along the z direction (resulting in a macroscopic dis- 



tortion of the crystal), whereas the first component of E u 
stretches tetrahedra A and squeezes tetrahedra B along 
the same axis (leaving the crystal dimensions unaltered 
to leading order). The resulting distortions of tetrahedra 
A and B can be written 

x A = (x g + x„)/\/2, x s = (x g - x„)/V2. 

The sum of elastic and magnetic energies, 

& = — J [lA ■ Xa + Ib • XflJ H 1 , 

is readily minimized with respect to the phonon variables: 

E ^ = 4k g (10) 

The minimized energy thus consists of two terms. The 
first term is diagonal in £4 and L3: 

-J' (^ 1 + fc- 1 )(/i + /l)/4. 

It puts tetrahedra of both types into one of the three 
collinear states (red, green, or blue), thus defining a 3- 
state Potts model. The cross term, 

introduces a coupling between the Potts states on the two 
sublattices. A softer even phonon (k g < k u ) yields the 
ground state of a ferromagnetic Potts model: all tetrahe- 
dra have the same primary color. A softer odd phonon 
(kg > k u ) produces a ground state with two different 
primary colors on the two sublattices. 

Translated back into spin language, the two ground 
states are shown in Figs. || and The latter, in fact, 
describes the Neel state observed in YM112 and MgV2C>4, 
compounds with spontaneous structural distortions. 

III. LANDAU THEORY 

Our simple model of classical spins on an elastic lat- 
tice of tetrahedra appears to be reasonably successful 
in explaining ground-state properties of some frustrated 
magnets. Can we also gain some understanding of phase 
transitions in these materials? 

To start with, we need to identify the relevant phases. 
At high temperatures, we have a symmetric paramag- 
netic state with no spin or bond order and no lattice dis- 
tortions. The ground state (T = 0) is a Neel phase with 
a distorted lattice. The two phases are distinguished, 
for example, by spin averages (Si) and by the disparities 
in bond lengths. In general, there may (and in certain 
cases will) exist an intermediate spin-Peierls phase. It 
is distinct from the Neel phase by the absence of spin 
order ((S,) =0). It is also different from the param- 
agnetic phase by the presence of lattice distortions and 
unequal spin correlations (Si-Sj) between various near- 
est neighbors, with a concomitant lowered symmetry. In 
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PIG. 7: Spontaneous dimerization of a spin chain. 

such cases, we expect two phase transitions: first, a high- 
temperature spin-Peicrls transition, which partially re- 
lieves frustration of spins, then, at a lower temperature, 
a transition into a Neel state. Such a scenario is permit- 
ted by symmetry and the frustration makes it easier to 
generate fluctuations that destablize the Neel state with- 
out destroying the bond order. 

In this section, we will discuss the spin-Peierls transi- 
tion using Landau theory. By analogy with spin chains, 
we will identify the relevant order parameter and discuss 
possible phase transitions in the framework of the Lan- 
dau theory. 

A. Dimerized chain 

To introduce method and notation, we start with a 
familiar example, namely the Landau theory of the spin- 
Peierls transition in antifcrromagnetic chains coupled 
to three-dimensional phonons. The paramagnetic and 
dimerized phases can be distinguished using the order 
parameter 

/ = (S2n+l'S2n — S2„'S2„-l). (11) 

It vanishes in the paramagnetic phase, since all bonds 
are equivalent. Spontaneous dimerization increases the 
probability of finding a singlet on half of the bonds, which 
leads to a nonzero value of /. Expansion of the free en- 
ergy (per spin) in powers of the order parameter contains 
even powers of / only: 

F(f, T) = F(0, T) + a(T)f 2 + c(T)/ 4 + . . . (12) 

Assuming that a becomes negative below T = T c , so that 
a(T) sa a(T — T c ), while c > and roughly constant, one 
obtains the standard scenario of a second-order phase 
transition: the minimum of the free energy shifts con- 
tinuously from / = above T c to / = ±yJa\T - T c |/2c 
below T c . 

The continuity of the transition depends crucially on 
the absence of a cubic term in the expansion (O) . With 
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chains, its absence is guaranteed by symmetry: states 
with / differing only by a sign are physically equivalent 
(Fig. 0), hence only even powers of / are allowed. 

Formally, the fate of an / 3 term is decided by its sym- 
metry properties. The symmetry group of an undistorted 
chain includes inversion on a site, which takes / i— > — /. 
Likewise, / 3 i— » — / 3 . Since, however, free energy must 
be invariant under all symmetry transformations, an / 3 
term is forbidden. 

In contrast, we will find that a cubic term is allowed 
in certain cases for the spin-Peierls order parameter on 
the pyrochlore lattice. In such cases, the spin-Peierls 
transition is expected to be discontinuous. 



B. Pyrochlore lattice: order parameter and broken 
symmetries 

What order parameters would characterize a spin- 
Peierls phase in a network of tetrahedra? The small- 
est unit of the lattice, a tetrahedron, contains 6 bond 
variables, so that there are 6 averages (Si • Sj) and 5 
differences between them, all of which could serve as or- 
der parameters. From the symmetry viewpoint, they can 
be divided into irreducible representations of the tetra- 
hedron group. One of them is the doublet f = (/i,/2), 
where 



h = 



h = 



<(Si + S 2 )-(S 3 + S 4 ) - 2Sx • S 2 - 2S 3 • S 4 ; 
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((Si - S 2 ) • (S 3 - S 4 )> 



(13) 



The other is a triplet - see Eq. (||). In the paramagnetic 
phase, both the doublet and triplet order parameters 
must vanish (all nearest-neighbor bonds have the same 
strength). In a spin-Peierls phase, either the doublet, or 
the triplet (or, potentially, both) will have nonzero expec- 
tation values. Energy considerations of Section |l| suggest 
that the driving force of this transition is the doublet. 

The two-component order parameter f can be the same 
for all tetrahedra, in which case only the rotational sym- 
metry of the lattice will be broken. Symmetry with re- 
spect to inversion on a site can also be violated if the 
order parameter f is not the same on tetrahedra of dif- 
ferent orientations. Lastly, translational symmetry of the 
lattice can also be broken if f varies among equivalent 
tetrahedra forming a commensurate wave. 



C. Pyrochlore lattice: q = phonons 

We restrict the analysis to situations when the trans- 
lational symmetry of the lattic e remains intact, as we 

In this case, any two 



did previously in Section [I B 



tetrahedra of the same orientation distort in the same 
way reducing the space group of the pyrochlore lattice 
to the octahedral point group Oh = I x (inversion / 



exchanges tetrahedra of different orientations). Despite 
this rather drastic simplification, we will see that we can 
account for the experimentally observed behaviour of a 
number of compounds, at least qualitatively. Phonons 
at other points in the Brillouin zone may drive a spin- 
Peierls transition as well, leading to ordered states with 
larger and often more complex unit cells. 

The order parameter has (potentially unequal) val- 
ues fyi and ts on tetrahedra of inequivalent orienta- 
tions. Their symmetric and antisymmetric combinations 
g = + fa and u = fA — ffl are irreducible doublets of 
the group Oh- In the paramagnetic phase, g = u = 0. 
In various spin-Peierls phases, one or both of these order 
parameters are nonzero. 

For classical spins, the domain of possible values of 
the order parameters f = (/i,/2) is the familiar tri- 
angle shown in Fig. ^. In view of the three-fold sym- 
metry (more precisely, permutation group S3), the two- 
dimensional vector f can be interpreted as a color. The 
extremal points represent the red, blue, and green states 
with collinear spins. In the paramagnetic (white) state 
f = 0, the color symmetry is manifest: the three primary 
colors are represented equally. In any spin-Peierls phase, 
the global color symmetry S3 is spontaneously broken. 



1. Landau free energy 

The Oh symmetry of the high-temperature phase al- 
lows the following terms in the Landau free energy: 



F(g, u) = a g g 2 + b g g 3 cos W g + c g g A 
+ a u u 2 + c u u A + d u u 6 cos 69 
+ b u u 2 g cos {26 u + 6 g ) + ... 



(14a) 
(14b) 
(14c) 



Here g = {g cosO g , gsm6 g ) with an analogous definition 
of u and 9 U . The first (second) line contains the lead- 
ing terms for the even (odd) distortion; the third line 
represents the lowest-order coupling between g and u. 
The constants a through e in this expression cannot be 
determined by symmetry considerations alone; when con- 
venient, one can try to determine their likely sign by tak- 
ing recourse to microscopic model Hamiltonians for the 
spin-lattice system. 

Omitted higher-order terms are assumed to be positive 
for stability. Landau theory is of course strictly to be 
applied only for small values of the order parameters. 
However, the shape of the range of the vector f (Fig ^) 
encodes some information on where the order parameter, 
once close to saturation, may point. 

This form of the free energy permits a number of dis- 
tinct ordered states. Generally, the symmetry of the lat- 
tice is reduced from cubic to tetragonal. In addition, the 
presence of an odd distortion (u ^ 0, or ^ fs) also 
breaks the symmetry of inversion through a site, exchang- 
ing tetrahedra A and B. Note that, whene ver a staggered 
distortion u is present, the coupling term (14c) generates 
a subdominant uniform distortion g of the crystal. 
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FIG. 8: Development of the order parameters and fs in the 
six scenarios of the spin-Peierls phase transition, f = in the 
paramagnetic (white) state, (a-b) A first-order transition is 
driven by the E g phonon and is discontinuous, = fs 7^ 0. 
(c-d) As a result of that transition, the E u phonon softens 
and triggers a subsequent second-order transition into a phase 
with f/t 7^ fs • (e-f ) A second-order transition is driven by the 
E u phonon. The E g order parameter is also induced, so that 
U + -fs- 



The phase transitions can be first or second order, de- 
pending on the mode driving the transition: the free en- 
ergy of the even mode g may have a cubic term ( |14aj ), 
which generally leads to a discontinuous jump . Th e odd 
mode u does not have its own cubi c te rm (14b), but 
is instead coupled nonlinear ly to g (14c). This differ- 
ence has the following physical origin. When an even 



distortion is present, all tetrahedra have, say, 4 strong 
and 2 weak bonds (a state of primary color). Chang- 
ing the sign of the order parameter g would give a state 
with 2 strong and 4 weak bonds on every tetrahedron 
(secondary color) , which need not have the same free en- 
ergy, F(g, 0) 7^ F(— g, 0). Hence a g 3 term is allowed. 
On the other hand, in a state with a pure odd distor- 
tion, tetrahedra A and B have opposite colors (e.g., A 
is blue and B is yellow) and switching the sign of u 
merely exchanges them (A is yellow and B is blue), so 
that F(0, u) = F(0, — u). Hence, there is no u 3 term. 



2. Spin-Peierls phases 

In the high-temperature paramagnetic phases, a g > 
and a u > 0, the minimum of the free energy lies at g = 
u = £4 = fs = 0. At low enough temperatures, one 
or both of these coefficients may become negative — see 
Eq. ([To|) . The nature of the resulting phase transition 
depends sensitively on the order in which a g and a u turn 
negative, as well as on the signs of the Taylor coefficients 
b u , b g , d u . Our results for the Heisenberg model (Section 
U) are compatible with the choice b u > 0, which we will 
assume in what follows. Below we describe six scenarios 
depicted in Fig. |L 



(a) In the simplest case, the even mode g becomes 
unstable, while the odd mode remains suppressed at all 
temperatures. The transition is disc ontin uous because 
of the cubic term in the free energy (14a). For b g > 0, 
minima of the free energy are at = 7r, ±7r/3. Asu = 0, 
distortions are the same on all tetrahedra, = fg = 
g/2. Thus each tetrahedron shows the same tetragonal 
distortion with 4 strong and 2 weak bonds (one of the 
three primary-color states). 

(b) Same as (a) but b g < 0. The minima are at 
6 g = 0, ±27r/3, secondary-color states with a tetragonal 
distortion, 2 strong and 4 weak bonds on all tetrahedra. 

(c) As the even order parameter g grows, it modifies 
the quadratic term of the odd mode through the non- 
linear coupling (14c). Once a u — b u g vanishes, scenarios 
(a) and (b) are modified: a second, continuous transi- 
tion occurs into a state where both g ^ and u / 0, so 
that fji 7^ fs. The directions of g and u are correlated: 
26 u + 9 g — it. For b g > (a), u is parallel to g; therefore, 
vectors and fs still point towards one of the corners of 
the triangle, but their length differ. Distortions of tetra- 
hedra A and B are remain tetragonal, but are unequal in 
strength. The symmetry of inversion is broken. 

(d) When the odd mode softens for b g < (b), u is 
perpendicular to g. The uniformly distorted state (e.g., 
yellow) becomes nonuniform: tetrahedra A acquire a red 
component, tetrahedra B acquire a green one. Distor- 
tions of individual tetrahedra are no longer purely tetrag- 
onal: there is an orthorombic component. Because the 
latter has a staggered nature, the lattice as a whole re- 
tains tetragonal symmetry. The symmetry of inversion 
is broken. 

Caveat. Because the high-temperature transitions in 
cases (c) and (d) are discontinuous, the intermediate 
phase (single primary or secondary color) may be skipped 
completely. In that case, instead of a succession of two 
transitions, there will be a single, discontinuous transi- 
tion directly into the final state with two different colors, 

(e) The transition can also be driven by the odd 
phonon, in which case it is expected to be continuous. 
The initial direction of the vector u is determined by 
the sign of the sixth-order anisotropy d u . For d u < 0, 
6 U = mr/3, where n is an integer; vectors and fs 
point in opposite directions, say, (a has a primary color 
(e.g., blue), while fs has a secondary one (in this case, 
yellow). The nonlinear coupling term (14c) generates a 
subdominant order parameter g = 0(u 2 ) parallel to u. 
This parasitic order parameter increases the length of 
the primary-color component (£4 becomes a deeper blue) 
and reduces that of the secondary-color component (fs 
is a pale yellow). The distortions are tetragonal but of 
oposite directions on tetrahedra A and B. 

(f) Odd phonon with d u > 0. The free energy ( p.4b| ) 
has a minimum for 8 U = (2n+l)7r/6. Initially, distortions 
of tetrahedra A and B are orthorombic, e.g. along f% (a 
red or green, respectively, with a touch of blue). The 
parasitic component g = 0{u 2 ) and perpendicular to 
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u, bends fA and fs towards the primary red and green 
directions and makes the individual distortions mostly 
tetragonal (along orthogonal axes in real space). Note 
that the final state is the same as in (d). 

D. Relation to 3-state Potts models 

As we have already mentioned, the symmetry of the 
bond variables f (permutation group S3) iayokes a simi- 
larity to the Potts model with q — 3 statesO with energy 

E=jJ2K*i, (15) 

(ij) 

S{ = f,2,3 being Potts states. Indeed, a similar two- 
component order parametep-has been introduced for the 
q = 3 Potts model by Ono.E.3 The pure Potts states cor- 
respond to collinear spin configurations. In the current 
context, the order parameter lives on tetrahedra, which 
form a three-dimensional diamond lattice. It is entirely 
plausible that the spin-Pcierls transitions described in 
this paper should be analogous to phase transitions in 
Potts models with short-range interactions. 

To this end, we can identify the simplest scenario 
[Fig. ||(a)] with the ferromagnetic Potts model. The 
latter is known to have a first-order transition in three 
dimensions, E3 which is consistent with our mean- field re- 
sult. 

Transitions shown in Figs. 0(e) and (f) have their 
analogs in the antiferromagnetic 3-state Potts model. Re- 
sults (mostly numerical) for latti ces , in d = 3 dimensions 
have been obtained fairly recentlyOEjO Banavar et alx3 
have studied the model on the simple cubic lattice and 
found at low temperatures an ordered state with a bro- 
ken sublattice symmetry (BSS). As the name suggests, 
the two sublattices are inequivalent: spins on one sub- 
lattice are primarily in one Potts state (say, blue), while 
the other sublattice is dominated by the remaining spin 
states (red + green) in equal proportions (yellow). More 
recently, Rosengren and Lapinskas predicted the exis- 
tence of another phase, with permutationally symmetric 
sublattices (PSS), where sublattices. A and B are domi- 
nated by two different Potts states,E3 e.g., red and green, 
respectively. Their Monte Carlo simulations suggest that 
the 3-state Potts antiferronLaffnct on the diamond lattice 
orders into the PSS phase.ll3 In both cases, the phase 
transition appears to be continuous, with critical prop- 
erties of the XY model in d = 3. 

IV. NEEL PHASES 

The spin-Peierls transition, whether in chains or in 
three-dimensional magnets, is driven by the desire of 
spins to reduce frustration. In the bond-ordered phase, 
exchange strength varies from bond to bond because of 
the distortion. Thus frustration is relieved and the clas- 
sical ground state becomes unique, modulo global spin 



rotations. In three dimensions, we can expect a spin- 
ordered state at zero temperature. As argued before, the 
transition into a Neel state need not coincide with the 
spin-Peierls transition. Therefore, generally there will be 
three separate phases: paramagnetic, spin-Peierls, and 
Neel. (In those cases when the spin-Peierls transition is 
discontinuous, the system may go directly into the Neel 
phase bypassing the spin-Peierls stage.) 



A. Neel orders 

Particulars of the Neel order on a distorted lattice ob- 
viously depends on the details of the distortion, which 
strengthens some bonds and weakens others. Because 
precise knowledge of spin interactions is rarely available 
(even for the undistorted state!), one can try an alterna- 
tive route, namely to include spin averages (Si) in the 
Landau theory developed above for a spin-Peierls phase. 

To keep technical details to a minimum, we will restrict 
the discussion to Neel states that do not break transla- 
tional symmetry of the crystal, i.e., spin averages (Sj) 
will be assumed to be identical for all tetrahedra of the 
same orientation. Put another way, (Sj) is the average 
spin on the i-th sublattice, i = 1 ... 4. Evidently, this 
parametrization adequately describes only a fraction of 
possible antiferormagnetic orders. For example, one of 
the Neel ground states obtained in our simple magnetoe- 
lastic model (Fig. ||) is already beyond its scope. More 
generally, even q = bond ordered states with different 
strengths on the pair of bonds of inequivalent tetrahedra 
related by inversion cannot be translated into a q = 
spin state. 



1. Undistorted lattice 

Let us construct the Landau free energy for spins on an 
undistorted lattice. Using the order parameters Sj = (Sj) 
one obtains the following expansion for the free energy: 

i=l \i=l / i=l 

(16) 

In an antiferromagnet, a\ > 0; for stability, we take 
bo > 0. When ao becomes negative, the minimum of 
the free energy shifts away from Sj = and the system 
will enter a Neel state. The free energy is minimized by 
any configuration of spin averages satisfying J2i=i s « = 0j 
the length of the averages is given by sf = ~a^/2b^. The 
Neel pattern is thus not unique, as expected for a frus- 
trated magnet. 

In addition to the quartic term shown in Eq. (|l9|), the 
free energy expansion may contain one more quartic in- 
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FIG. 9: Neel order in the presence of a tetragonal distortion, 
as given by Landau theory. The abscissa is the distortion 
amplitude (a2/ao)/i; the ordinate is the biquadratic coupling 
61/60 • The four Neel phases are described in the text. 
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FIG. 10: A q = Neel state with orthogonal spins. The 
distortion is created by the same E g phonon as in Fig. ^ but 
with the opposite sign (blue bonds are enhanced). 



variant, 



i>j 



(17) 



This term, in fact, will break the degeneracy of the Neel 
states. In the case 61 < 0, the Neel phase has collinear 
spins (any one of the three collinear states in Fig. ^. 
Note that these states also break the bond symmetry.) 
When b\ > 0, the spin averages point at equal angles of 
arccos(— 1/3) ~ 109° to one another. 



2. Distorted lattice 

A lattice distortion, however small, breaks the cu- 
bic symmetry, so that additional, less symmetric terms 
will appear in the free energy. For a distortion that, 
symmetry- wise, belongs to the irreducible doublet E, the 
lowest-order perturbation will be of the form 



a 2 (f -0) = as (/i </>i + hfc), 



(18) 



where f is the familiar spin-Peierls order parameter de- 
scribing the distortion. The spin part 0({si}) should 
therefore also be a doublet of the same symmetry: 

01 = [(si + s 2 )-(s 3 + s 4 ) - 2si • s 2 - 2s 3 • s 4 ]/VT2, 

2 = (si - s 2 ) ■ (s 3 - s 4 )/2. (19) 

Apart from global rotations of all 4 spins, these two vari- 
ables uniquely determine the relative directions of the 4 
vectors Sj satisfying the constraint Y*t=i Si = ®- ^ * s con ~ 
venient to separate the direction and length of the spin 



averages in the free energy: = s * 
then takes the following form: 



0. The free energy 



F(s, 0) = a s 2 + 02 (f -0)s 2 + 60s 4 + 6i(0-0)s 4 . (20) 

The last term is simply Eq. fll7|); we have also dropped 
the term (X^=i s «) 2 assuming that the minimization is 
done over antiferromagnetic states. 

Minimization can now be done separately over the 
length s and direction variables 0. The minimization 
with respect to s at a fixed is straightforward giving a 
minimum 



F(0) = inf F{s,4>) 



[ao + a 2 (f-0)] 2 
2[6 o + 6i(0-0)]' 



(21) 



The minimization with respect to is done over the tri- 
angular domain shown in Fig. ||. The outcome is decided 
by a competition between the biquadratic exchange ( |T7j ) 
and the coupling to the spin-Peierls order (|l8|). Gener- 
ally, a distortion f pulls the vector in the same direction 
(if ao < and a 2 < 0), whereas the biquadratic coupling 
attempts to minimize (61 > 0) or maximize (b\ < 0) its 
length. 

Fig. |^ shows the appropriate phase diagram of the anti- 
ferromagnetic ordering for the case of a uniform tetrago- 
nal distortion. When the influence of the distortion dom- 
inates, we find two Neel phases. For a distortion that 
produces four strong bonds per tetrahedron, the spins 
are collinear (e.g., si = s 2 = —S3 = — s 4 ); a distor- 
tion of the opposite sense (four weak bonds) stabilizes a 
coplanar state with two orthogonal pairs of spins (e.g., 
Si = — s 2 , S3 = — s 4 , while Si-s 3 = 0, Fig. Il0|). In an 
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intermediate state, where the biquadratic coupling dom- 
inates, the spins are no longer coplanar. For a positive 
biquadratic coupling 61, they gradually interpolate be- 
tween the collinear and orthogonal orientations as the 
strength of the distortion varies. In the shaded region 
(bi < 0), a uniform tetragragonal phase becomes unsta- 
ble: the distortion acquires an orthorombic component 
at the onset of the Neel order; in addition, spin averages 
Si break the translational symmetry of the crystal, as in 
Fig.|. 



V. CONCLUSION 

In this paper, we have presented a theory of the spin- 
Peierls transition in a frustrated magnetic system, the 
Heisenberg antiferromagnet on the "pyrochlore" lattice. 
Several aspects distinguish this effect from its counter- 
part in spin chains (e.g., CuGeOa). (1) The magnetic 
system is manifestly three-dimensional, initially possess- 
ing a cubic symmetry (2) The effect is classical: the 
quantum mechanics of spins plays no significant role. (3) 
The order parameter (disparity of spin correlations) has 
rather nontrivial symmetry properties: its components 
form an irreducible doublet of the tetrahedral symmetry 
group. 

What drives this spin-Peierls transition? From the 
kinematical viewpoint, the color variables (^) can be 
viewed as coordinates in the manifold of ground states. 
The onset of a bond order thus lifts the large accidental 
degeneracy of the ground state. The dynamical reason 
for the transition is the Jahn- Teller effect occuring in the 
building blocks of the "pyrochlore" lattice, tetrahedra of 
magnetic ions. Unlike in many cubic spinels, where the 
Jahn- Teller distortion is caused by an orbital degener- 
acy, here it is driven by the spin degrees of freedom. An 
isolated tetrahedron (whether with quantum or classical 
spins) undergoes a tetragonal distortion along one of the 
three orthogonal symmetry axes, producing four strong 
and two weak bonds — or vice versa, depending on the 
sense of the distortion. (The residual threefold degen- 
eracy of the spin-and-phonon ground state relates this 
problem to the 3-state Potts model.) 

A study of the spin-driven distortion of an isolated 
tetrahedron provides clues to the simplest description 
of the spin-Peierls effect on the lattice of corner-sharing 
tetrahedra, in particular the form of the order parameter. 
Depending on the parity of the phonon responsible for the 
transition, Landau theory predicts a first or second-order 
spin-Peierls transition. It is gratifying to see that these 
predictions are consistent with known properties of the 
3-state Potts model on the diamond lattice, onto which 
our problem can be mapped. 

Although a distortion of the lattice reduces the geomet- 
ric frustration, there is no reason to expect that it induces 
a Neel order immediately at T c . With the exception of 
strongly discontinuous spin-Peierls transitions, spin or- 
der will will generally set in at a lower temperature than 



bond order. A transition between a spin-Peierls and Neel 
phases has been analyzed above, again in the framework 
of Landau theory, which predicts collinear, coplanar, and 
more general antifcrromagnetic orders at the lowest tem- 
peratures. 

Are there experimental realizations of the bond- 
ordering transition described in this work? The spinel 
compound ZnCr204, in which Cr 3+ ions (S = 3/2) form 
the tetrahedral network, appears to be a good candidate 
for a Heisenberg antiferromagnet on the "pyrochlore" lat- 
tice. Indeed, it exhibits a magnetic transition. a ^ T c = 12 
K accompanied by a structural distortion. Ej'Ej Neel order 
is observed at T < T c , which is not surprising given that 
the transition is discontinuous. Another spinel, MgV204, 
shows a sequence of two transitions^ a structural one 
occurs at T C 2 = 65 K, Neel order appears at T c % = 42 
K. However, because of an orbital degeneracy (the outer 
electrons in V 3+ are 3s 2 p 6 d 2 ), the upper transition may 
well be triggered by the ordinary Jahn- Teller effect com- 
mon to spinels. Therefore, we cannot positively identify 
it as a spin-Peierls transition. Nevertheless, the Neel or- 
der in this compound (and also in YMn2) agrees with 
the prediction of a simple magnetoelastic model given 
here (Fig. ^). Lastly, a recently discovered second-order 
structural phase transitionEa'Ell in the metallic pyrochlore 
Cd2Re207 at T c = 194 K may turn out to be a closely 
related Peierls transition, whose symmetry properties are 
identical. 

We have presented a bare-bones theoretical description 
of the spin-Peierls transition in a "pyrochlore" antiferro- 
magnet. Only the simplest ordering patterns have been 
discussed, namely those that do not break the transla- 
tional symmetry of the crystal. An obvious extension of 
this work would be to include bond and spin orders at 
nonzero commensurate wave vectors. E.g., a phonon with 
q = (it, 7r, 7r) appears to be responsible for the distortion 
in ZnCi'204. 
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APPENDIX A: A MODEL WITH A COPLANAR 
GROUND STATE 

In addition to pairwise spin exchange, which gives rise 
to the Heisenberg interaction Sj ■ Sj, spins can be in- 
volved in longer exchange cycles, such as 123 1— > 231 and 
312. The cyclic exchange of 3 spins induces the same 
pairwise Heisenberg interaction, which has already been 
considered. The next nontrivial contribution comes from 
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FIG. 11: (a-c) Cyclic exchange of 4 spins, (d) A ground state 
with two orthogonal pairs of (classical) spins. 



4-spin cyclic exchange 

P1234 = (S 1 -S 3 )(S 2 .S4) + (S 1 -S 4 )(S 2 .S 3 ) 
- (S 1 -S 2 )(S 3 -S 4 ), 

which — for S = 1/2 — moves spin states clockwise or 
counterclockwise around the loop 1234 (Fig. |ll|). For 
localized spins, this interaction is weaker than pairwise 
exchange and can be considered as a perturbation. ilev- 
ertheless, its signature has apparently been detected^ in 
the spin- wave spectrum of La 2 CuC>4. 

A tetrahedron has three loops of length four: 1234, 
1324, and 1243. Quantities -P12347 P1324, and -Pi 2 43 con- 
tain the trivial representation A\ and the doublet E. 

The sum Pi 23 4 + Pi 3 24 + -Pi 24 3, being invariant under 
all symmetry operations, can enter an expression for the 
energy on its own. In the subspace of the ground states 
of the Heisenberg Hamiltonian, this perturbation can be 
written as a biquadratic term 



i>j 



S 3 ) 2 = 



3J4/ 2 



const. 



The outcome depends on the sign of the coupling con- 
stant J4. If the constant is negative, it prefers the largest 
magnitude of f and selects the three collinear spin states 
(primary colors in Fig. ^J). A positive J4 suppresses f, 
selecting the white state with all bonds equivalent; spin 
averages make angles of arccos(— 1/3) s» 109° with each 
other. 

The two differences between Pi2 3 4, Pi 3 24, and P i2 4 3 
form an irreducible doublet E. They can therefore couple 
to the phonon doublet of the same symmetry adding this 
magnetoelastic term to the energy: 



Ja 



1 ( -Pl324 + -Pl243 — 2P i23 4 -Pl324 — Pl243 



V6 



-Xi 



^/2 



-;r 2 



where J4 describes variation of the cyclic exchanges un- 
der a tetragonal or orthorombic distortion. After adding 
an elastic term kx 2 /2 and minimizing with respect to the 
phonon variables x, one obtains the following contribu- 
tion to the energy: 



2J' A 



3k 



L{ [( Sl .S 2 ) 2 -(S 2 .S 3 ) 2 ] 2 

[(S 2 .S 3 ) 2 -(S 3 -S 1 ) 2 ] 2 
[(S 3 -S 1 ) 2 -(S 1 -S 2 ) 2 ] 2 }. 



(Al) 



The energy is lowered by the greatest amount in coplanar 
states with spins making angles of 90° and 180° with one 
another. These states have two pairs of frustrated bonds 
and are marked in secondary colors in Fig. ||: e.g. a state 
with blue and red frustrated bonds is magenta. 

As the two physical forces described above — the four- 
spin exchange and its coupling to the phonons — favor 
different ground states, the outcome is decided by their 
relative strengths. In particular, the ground state has 
orthogonal spins [Fig. [y](d)] when the 4-spin exchange 
is very sensitive to atomic displacements (and therefore 
the distortion is large). The reason for this effect can 
be understood as follows. A strong tetragonal distortion 
enhances the 4-spin exchange around one loop and sup- 
presses the same around the remaining two. When the 
latter are completely switched off, we are dealing with a 
square, and for a square the 4-spin cyclic exchange pro- 
duces a ground state with orthogonal spins .EEl 



APPENDIX B: SINGLE TETRAHEDRON: 
QUANTUM SPINS 

Quantum spins coupled to classical distortions of the 
tetrahedron have the following Hamiltonian in the sub- 
space of ground states: 



H = -J'i-x + kx 2 /2, 



(Bl) 



which is formally the same as the classical energy (|7|). 
For a fixed "direction" of the distortion 

n = tc/x — (cos a, sin a), 

the operator f • n = j\ cos a + / 2 sin a has 2S+1 eigenval- 
ues A CT , <j — ... 25* in the ground-state manifold. The 



Hamiltonian (Bl) has the following energy levels: 

E a (x, n) = -J'\ a x + kx 2 /2. 

Minimization with respect to the magnitude of the dis- 
tortion x gives the following result: 

E a (n) = inf EJx,ii) = -J' 2 A 2 /2fc, 

X 

The energy is lowest in a spin state with the largest (in 
absolute terms) eigenvalue A CT of the operator f • n. The 
final step is minimization with respect to the direction n. 

Let us illustrate the mini mizat ion procedure using the 
classical problem of Section II A as an example. To this 
end, we show classically allowed values of the vector 
An = (f • n) n as shaded areas in Fig. |l2| (the values of 
f fill the interior of the regular triangle shown in dashed 
lines). The largest magnitude |An| = |A| is found in the 
directions a = 0,±27r/3, which correspond to the three 
collinear states. In the quantum case, A CT has a discrete 
spectrum, therefore allowed values of Ao-n will show as 
lines on the same graph. 
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FIG. 12: Solid lines: eigenvalues \ a of the operator f-n = 
fi cos a + fi sin a in polar coordinates ( X a , a) for several spin 
lengths S. Shaded area: the classical result for S — > c« (prop- 
erly rescaled). 

The main subtlety of the quantum problem is the non- 
commutativity of the bond operators f\ and / 2 : 

[/i,/ 2 ] = -2V3i X) (B2) 



where \ is the operator of chirality 

X = S 1 .(S 2 xS 3 ) = -S 1 .(S 2 xS 4 ) = ... 

Therefore eigenvalues of f • n cannot be constructed from 
those of fi and / 2 , but rather must be determined for 
every direction fi. 



1. Matrix elements of the operator f 

In the (25 + l)-dimensional subspace of singlet ground 
states we choose the basis {|cr)}, a = ... 25 being the 
total spin of pair 12. (The spin of pair 34 must be the 
same in order to form a total spin of 0.) The operator f\ 
is also diagonal in this basis because 

fx = [Si 2 • S 34 - S? 2 - S^ 4 + 45(5 + 1)]/Vl2. 

Its eigenvalues are 

/i = [45(5+1) -3<r(<7 + l)]/Vl2. (B3) 

The operator / 2 is off-diagonal. To compute its matrix 
elements, write out an expression for the singlet \a): 

\J la + 1 — 

v fj,=—a 

Here \a, n)\2 is the state of pair 12 with total spin a and 
its projection fj, onto a chosen axis. Then, by definition, 



(7 <7 



(a'\h\a) = - = J2 E {-iy + °'-^' (a',VI(S3-S 4 )|a,- / ,)34-(a',M'|(S 1 -S 2 )|a, / ,) 12 

2 V( 2(T + -LJ(2cr+ Ij l= _ a , a 



<7' U 



2 v .2 ( 7' + l)(2a+l) „£,,£ (CT,Ai 



Si - S 2 )|<r , (j, )i2 • {a , ii | (Si - S 2 )|cr,/z)i 2 



(B4) 



— — a' /U=— u 



r 



In simplifying this expression, we have replaced indices 3 
and 4 with 2 and 1 because the matrix elements involve 
one pair at a-tjme. Also, we have used the properties of 
time reversalc3 to simplify the first matrix element in the 
summand. 



The last line of Eq. (B4) contains matrix elements of 
Si — S 2 , which is (a) a vector; (b) antisymmetric in 1^2. 
These lead to a selection rule: 



Completeness of the basis {|cr',/i')} allows us to further 
simplify the right-hand side at the expense of producing a 
set of 25+ 1 coupled equations. This is done by summing 
over a' with appropriate weights: 



(o-'l^k) = unless a' = a ± 1. (B5) 
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2S 

E 

cr'=0 



2a' 



2a 



Y^'l/2k) = E (^KSi - S 2) 2 k: M>i2 = 25(5 + 1) - + 1). 



Solving these gives the matrix elements of f^. 

The operator f has the following nonzero matrix ele- 
ments: 

45(5+1) -3^ + 1) 
(ol/iF) = t= , (B6) 



12 

(a - ll/ak) = W/ak - 1) = "^y// 21 . (B7) 



2. Eigenvalues of the operator fn 



For 5 = 1/2, we have 



„ „ \/3 / cos a 
f n = — — 

2 V sma 



sma 
cos a 



(B8) 



The eigenvalues A = ±^3/2 are independent of the di- 
rection n. Thus there is no preferred direction for the dis- 
tortion, unless one introduces some additional, nonlinear 
couplings. n This fact was noted previously by Yamashita 
and Ueda.Q 
For 5 = 1, 



4 cos a 4 sin a 
f ■ n = — ^ | 4 sin a cos a \/5 sin a 
y/E sin a —5 cos a 



(B9) 



Its eigenvalues are given by the equation 



X A - 7A 



20 cos 3a 
3^3 



0. 



The largest eigenvalue is attained when a = n or ±7r/3, 
i.e. when n points towards one of the corners of the 
triangle (the three tetragonal distortions). For 9 = tt, 
the total spin of bond 12 is a = 2 (the same goes for 
bond 34). For the other choices of a, the highest spin is 
found on other pairs of opposite bonds. 

We have checked higher spin values 5 < 10 and always 
found that the lowest energy is obtained for a tetrago- 
nal distortion, when two opposite bonds have the highest 
spin 25* and are thus most strongly frustrated. One can 
use perturbation theory to show that a — n is at least a 
local maximum of |A2s| for any 5 > 1/2: 



A 2 s(a) 2S-1 
A 2S (0) 45-1 



a- 2 + 0(a 4 ). 



As the value of S increases, eigenvalues X a fill the classi- 
cal region (three overlapping shaded circles in Fig. 12). 
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